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1. Introduction 

One of the most crucial features of QCD is the dynamical generation of a quark mass and 
the related phenomenon of chiral symmetry breaking (CSB). The realization of these characteristic 
QCD phenomena is intimately connected with the non-perturbative nature of the infrared sector 
of the theory; therefore, their study in the continuum is best studied within the inherently non- 
perturbative formalism based on the Schwinger-Dyson Equations (SDE) [|l|, §]. 

It is well known that the non-linear integral equation governing the behavior of quark prop- 
agator (gap equation) displays a very rich structure. In particular, the existence or not of non- 
trivial solutions for this equation depends crucially on the strength of its kernel, which is largely 
determined by the non-perturbative behavior of the fully dressed gluon propagator and the quark- 
gluon vertex [[I]]. Actually, the quark-gluon vertex plays an absolutely essential role, introducing 
to the gap equation a dependence on the ghost dressing function [||] and the quark-ghost scattering 
amplitude |[T]]; as we will see, this dependence turns out to be numerically crucial for obtaining 
phenomenologically acceptable quark masses. 

In this talk, we will show how it is possible to obtain physically relevant solutions for the quark 
dynamical mass using the infrared finite ingredients obtained from the lattice, i.e., a non-vanishing 
gluon propagator and a non-enhanced ghost sector. To do that, we build a truncation scheme for the 
quark SDE, supplemented with three non-perturbative ingredients: (i) the gluon propagator, (ii) the 
ghost dressing function obtained from large- volume lattice simulations [@, |5]], and (iii) the "one- 
loop dressed" approximate version of the scalar form factor of the quark-ghost scattering kernel [|l|]. 
In addition, the results for the fermion masses in the adjoint representation are briefly discussed. 

2. The gap equation 

Let us start by defining some basic quantities that are important for the analysis of the CSB. 
The quantity that will be on the focus of our attention is the full quark propagator; in Minkowski 
space, its inverse has the general form given by 

S-\p)=i>-m- L{p) =A{p')i> - B{p^) , (2. 1) 

where m is the bare current quark mass, and the quark self-energy. Notice that the self- 
energy can be decomposed in terms of a Dirac vector component, A{p^), and a scalar com- 
ponent, B{p^), which allow us to define the dynamical quark mass function as being the ratio 
^/{p^)=B{p^)/A{p^). 

As we are interested in generating the quark mass exclusively through dynamical effects, we 
consider the case without explicit CSB, where the bare mass m = 0. Then, it is easy to see that 
the quark propagator will develop a dynamical quark mass only if the scalar component, B{p^), is 
different from zero [^. 

The diagrammatic representation of S^^{p) is shown in Fig.|T| using the same convention of 
momenta flow as indicated in the figure, the gap equation can be written as 

S-\p) =^-Qg^ Jrfs{k)ry{-p,k,q)A^^{q), (2-2) 
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where q = p — k, Ji^ = ^■^^{In)^'' J d'^k, with d = A — e the dimension of space-time. Cr is the 
Casimir eigenvalue of the given fermion representation, where we set r = F for the fundamental, 
and r= A for the adjoint. Specifically, for the gauge group SU (3), we have Ca = 3 and Cp = 4/3. 



q=p-k 




Figure 1: Diagrammatic representation of the quark SDE (gap equation). 



According to Eq. (1.2), the quark propagator, S{p), is dynamically determined in terms of 



an integral equation involving itself, the full gluon propagator, to be denoted by A^^(^), the full 
fermion-gluon vertex rv{—p,k,q), and its tree level counterpart pji'^ = y^. 
In the Landau gauge the full gluon propagator, A^y{q), has the form 

A^^(^) = - 



r 



A{q^), (2.3) 



where the non-perturbative behavior of the scalar factor A{q^) has been studied in great detail in 
the continuum ||], and in the lattice simulations ^. 

2.1 The full quark-gluon vertex 

In principle, the fully-dressed quark-gluon vertex, Fy, is determined from its own SDE, which 
contains a number of (unknown) multiparticle kernels, characteristic of the so-called "skeleton 
expansion". Dealing with such an equation is technically very difficult; therefore, the standard way 
to obtain information about this vertex is to use a gauge-technique Ansatz for it [|l^]. The general 
idea of this method is to express the longitudinal part of a vertex in terms of the various quantities 
appearing in the fundamental Ward identity or Slavnov-Taylor identity (STI) that it satisfies. 

In the case of the quark-gluon vertex, the general Lorentz decomposition for the longitudinal 
part of Fy involves four different form factors. All these form factors contain an explicit depen- 
dence to both the ghost dressing function and the "quark-ghost scattering kernel" which, in turn, 
also has a rich tensorial structure. Specifically, the vertex r^{pi,p2,P3) satisfies the following 
STI[|] 

P3'^^{PUP2,P3)=F{p3)[S^H-pi)H{pi,P2,P3)-H{P2,PuP3)S^Hp2)], (2.4) 

where H{p\,p2,P3) is the quark-ghost scattering kernel represented in the Fig. ^ and F{p3) is the 
ghost dressing function which is related to the full ghost propagator by D{p3) = iF{p3)/ p^. 

The kernel H{pi,p2,P3) and the "conjugate" H{p2,pi,P3) have the following Lorentz decom- 



position [11] (note the change p\_ ^ P2 in the arguments of the latter) 



H{pi,p2,P3) = XoI + Xi^i+X2h+X30^vPiP2, 

H{P2,PuP3) = XoI-X2j^i -Xih+^i^nyPiP2, (2-5) 
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Figure 2: The fermion-ghost scattering kernel H{pi,p2,p3)- 



where the form factors X,- are functions of the momenta, Xi=Xi{pi,p2,P3), and we use the notation 
Xi{p,r,q) =Xi{r,p,q) and a^v = ^[yti,Yv]- 

The most general Lorentz decomposition for the longitudinal part of the vertex T^[pi,p2,p-i) 
can be written as [ITTll 



^ii{p\,P2,P?,) =Li7^ +^2(^1 -h){Pi -P2)ii+L^{p\ -P2)n+L40^v{pi -piY , (2.6) 

where L, are the form factors, whose dependence on the momenta has been suppressed, in order 
to keep a compact notation, i.e., Li = Li{p\,p2,P3). Notice that the tree level expression for pj^' is 
recovered setting L\ = \ and L2 = L3 = L4 = 0; then, rj^' = 7^. 
Contracting Eq. ( [2.6[ ) with p^, we have 

P3^fi{Pi,P2,P3) = {pl-p])L^I+[{pl-p'i)L2-Li]^i - [{pl-p\)L2 + L\\^2-2Ud^vp^^pl . 

(2.7) 

In addition, substituting the standard decomposition of S^^{p), given in Eq. (2.1), and the 
expression of H{pi,p2,P3) given by Eq. (2^) into Eq. {2A), we find that the rhs of Eq. (2J_) can 
be also expressed in terms of the functions A, B and X,'s. Then, it is relatively straightforward to 
demonstrate that the L, 's may be expressed as [p 

F{P3) 



+ 



L2 



+ 



2 
2 

FjPi,) 
2{pj-pj 

F{P^) 
2{pI-p\ 

F{P3) 

2 2 

Pi -pi 



A{p,)[Xo + (pj- PvP2)X3]+A{p2)\Xo + (pj- PVP2)X3 

B{pi){Xi +X2) +B{p2){Xi +r2)| ; 

A{pi)[{p\+pi ■p2)X3 - Xo] -A{p2) [{pI +PvP2)X3- Xo] 
B{pi){X2-Xi)+Bip2){Xi-X2)y, 
A(pi) {PIX1+PVP2X2) -A{p2) {plXi+pyp2X2) +B{pi)Xo-B{p2)Xo \ ; 



L4 



F{P3) 



A{p,)X2 -A{p2)X2+B{pi)X3-B{p2)X3 



(2.8) 



It is interesting to notice that setting in Eq. ( p.8| ) Xq = Xq = 1 and X,- = X; = 0, for / > 1 , and 
F{p3) = 1, we obtain the following expressions 



A{pi)+A{p2) 



A{p,)-A{p2) 
2{P\-Pl) 



L3 



B{p,)-B{p2) 



2 2 
Pi -Pi 



L4 = 



(2.9) 
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which give rise to the so-called Ball-Chiu (BC) vertex [12], 

A(pi)+A(/J2) n 



2 

-PiY 



-p\ 



+ [B{pi)-B{p2)] 



(2.10) 



which is widely employed in the literature for studies of CSB 

In this work, we go beyond the Abelian version of the quark-gluon vertex expressed by 



Eq. ( |2.10| ), and we consider the case where only the scalar component of the quark-ghost scat- 
tering kernel is non- vanishing i.e. Xq ^ while Xi = Xj = 0, for / < 1. In this limit, the expressions 



of Eq. (2.8) reduce to 



Li=f(p3)^0(P3) 
L2= F{p3)Xq{p3) 



A{pi)+A{p2) 



Aip,)-Aip2) 

M-pD 



L,=F{p3)Xoip,) 

u = o. 



B{p,)-B{p2) 

2 2 
P\-P2 



(2.11) 



According the above expression, the form factors L,'s display an explicit dependence on the 
product F {pt,)Xq{pt,) which contains information about the infrared behavior of the ghost propa- 
gator ||l]]. Therefore, the information about the ghost sector enters into the gap equation Eq. {1.2) 
through the quark-gluon vertex of Eq. (p^. 

Since the transverse part of the quark-gluon vertex is not constrained by the STI given by 
Eq. ( |2.4[ ), it is possible to add a transverse part to the vertex without violating the this important 
identity. The so-called Curtis and Pennington (CP) vertex [^], to be denoted by F^p, modifies Fj 
by an identically conserved term [0], such that 



BC 



(/J 1 , ;?2 , /?3 ) = r(^c (z' 1 ' /'z , P3 ) + (/? 1 , /'2 , ;?3 ) 



(2.12) 



with 



j^ill N y^{pl-p\)-{P\-P2Y{^l+h) ^ M ,0 1^, 

^t{P\,P2,P3) = „ ^ [A(P2) -A(pi)] , (2.13) 



where 



1 



ld{pi,p2) 



d{puP2) = 2 ] 2 { iP2-piY + 
P1+P2 



B^P2) , B^{p,) 



A^P2) A\p, 



(2.14) 



In analogy to what happens in the case of the BC vertex, the ghost effects due to F{p3) and 
the quark-ghost scattering kernel Xq\pj,) will be incorporated into the CP vertex through a simple 
multiplication of its tensorial structure by the factor F{p-i)XQ\pT,). Therefore, within our ap- 
proximation, the "ghost-improved" versions of the vertex, to be denoted by ^Yic{Pi^P2iP?,) and 
Tqy{p\,P2,P?,), respectively, can be written as 



^c{P\^P2-,P?.) = F{p3)X^'^{p3)^BciPUP2,P3 



^CpiPl^P2,P3) = F{p3)X^^\p3)r^p(Pi^P2^Pi)- 



(2.15) 
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Notice that the above equation in the limit of xl^^ {q) = 1 reduces to the vertex employed in Ref. [^]. 

Substituting the vertex T^c{PItP2,P^) given by Eqs. ( ^ ) and ( |2.1l[ ) into the gap equation 
(^^), and defining p\ = —p, p2 = k, and pj = q, we arrive at the following coupled system for 
A{p^) andB(/72) 



je^){p-k) 



A^{B)B + B^{B) 



Bip^ 



jeoip-k) 



(2.16) 
(2.17) 



where the kernel J^o(^) corresponds to the part that is not altered by the tensorial structure of the 
quark-gluon vertex, namely 



jeo{q)=mF{q)X^^\q), 
while the parts that are affected, je^^{k,p) and Jf^^{k,p), are given by [ |T4| ] 



(2.18) 



■y(fr{k,p) 



2p^ 

-A{k^)AA{k^,p^) 



'k^-^^ + 2'^hip,k) 
P P 



(2.19) 



^^''{k,p) = Uik^Mik^) +A{p^)]+2[B{k^)AA{k\p^) -Aik^)ABik\p')] h{p,k) , 



where 



and 



h{p,k) 



[k^p- - (k-p)- 



(2.20) 



(2.21) 



^2 _ p2 ' ' 1^2 _ p2 

Similarly, the effect of the vertex rQp{pi,p2,P3) is to replace the kernels J(f^^{k,p) and 
J(fg^{k,p), appearing in Eq. ( |2.16D and ( [2.17[ ), by J^^^{k,p) and Jf^^{k,p), respectively, where 



3k-p 



jf^'{k,p) = j(r^"'^{k,p) + ^-^A{k')hA{k\p', ^^^^^ , 

J^s^'ik^p) = .y(fi^{k,p) + ^Bik^)AA{k\p^) P^^' . 



J2\2 



(2.22) 



3. The non-perturbative ingredients from the lattice 

As we have seen in Eq. 2.18 , the gap equation given in Eqs. ( 2.16[ ) and ( |2.17| ) depends on the 
nonperturbative form of the three basic Green's functions, namely A(^), F{q), and Xo(^). There- 
fore, in order to proceed with the analysis, we use for A(^) and F(^) the recent lattice data obtained 
by [Q], and shown in Fig. |3[ 

We clearly see that both lattice results for A(^) and F{q) are infrared finite. This characteristic 
feature has been long associated with a purely non-perturbative effect, namely the dynamical gen- 
eration of an effective (momentum-dependent) gluon mass [ p^j , ^, |l^], whose presence saturates 
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Figure 3: Lattice results for the gluon propagator, A(^), and ghost dressing, F{q) renormalized at 
^ =4.3GeV. 

the gluon propagator in the IR. In addition, the appearance of the gluon mass is also responsible 
for the infrared finiteness of the ghost dressing function, F{q^) [0, [T7|], which is shown on the right 
panel of Fig. |3[ 

The only ingredient missing at this point is the form factor XQ\q). Unfortunately, as far we 
know, there are no lattice data available for Xo(^). We will therefore proceed with our analysis by 
deriving an estimate for Xq, based on its "one-loop dressed" approximation [|]]. 




Figure 4: Diagrammatic representation of the quark-ghost scattering kernel, H{~q/2, —q/2,q), at one-loop. 

Specifically, we will consider the scalar contribution of the diagram represented in Fig. ^ and 
can be written as 

X'^\p,p,q) = \-i(^^'^ h^\k)D{k - p) Gy Tr{S{k + q)r^ } . (3. 1) 

where the full gluon-ghost vertex is denoted by = d"^G^. 

Using the tree-level expression for Gy = {k — p)y, and the kinetic configuration where pi = p2 = p, 
and p = —q/2, we notice that the above equation simplifies considerably (for more details see [|l]|), 
and we arrive at 

xil\q) = l + ^CAgVl 



k^q^ 



A{k)F{k) 



F{k + q) 
ik + qf 



(3.2) 
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A non-perturbative estimate for Xq^ [q] can be obtained by substituting into Eq. (3^) the lattice 
data for A(g) and F{q), given in Fig. |3[ The numerical result for X^"^ (q) is shown in the Fig. 



1 







Form factor X'5(q) 

ffl Numerical result 
Fit 



























q'[GeV'] 



Figure 5: Numerical result for the form factor Xq ' (q) of the quark-ghost scattering kernel given by Eq. (3.2) 
when a (m^) =0.295. 

As we can see Xq^ (q) shows a maximum located in the intermediate momentum region (around 
450 MeV), while in the UV and IR regions xj^^ (q) ^ I. Although this peak is not very pronounced, 
it occurs in a region where the kernel of the gap equation is extremely sensitive. We have noticed 
that minor changes in the intermediate region result in considerable changes in the value of the 
dynamical mass generated as discussed in detail in Ref. 



4. Numerical results 



With all ingredients in hand, we are now in position to solve the system formed by Eqs. (2.16) 
and ( 2. 17| ). Substituting A{q-^), F{q^), andXQ''(^) into Eqs. ( 2.16 ) and ( 2.17 ), with the modification 
Z^^J^ B{k,p) — > J^ B{k,p)F{p^), to enforce the correct renormalization group behavior of the 
dynamical mass (see discussion in we determine numerically the unknown functions A{p^) 
and B{p^) for the fundamental representation. The result for the quark wave function A^^{p^) is 
shown in the left panel of Fig. ^, while the right panel shows the dynamical quark mass ^{p^). 

Notice that Fig. ^ shows the result for both vertices T^c{pi-iP2-,Pt) (red dashed curves) and 
^c?{P\iP2,P?,) (black continuous line). One clearly sees that ^(jp-) freezes out and acquires a 
finite value in the IR, (^(0) = 294 MeV in the case of the BC vertex, and ^(0) = 307 MeV for 
the CP vertex). Moreover, in the UV it shows the expected perturbative behavior represented by 
the blue dashed curve. 

Once the behavior of the dynamical quark mass has been determined, we have computed 
the pion decay constant and the quark condensate. For the pion decay constant we obtained 
fn = 80.6 MeV while for the condensate l^qq) (IGeV^) = (217MeV)^ which are in good agree- 
ment with phenomenological results [jl]]). 

Next, we solve the system of Eqs. ( [2.16| ) and (2. 17) in the adjoint representation i.e. Q = Ca = 3. 
When we switch from the fundamental to the adjoint representation, the overall effect in the gap 
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Figure 6: Left panel: The quark wave function A ' [p^) in the fundamental representation for both vertices 
Right panel: The dynamical quark mass ^{p^). 



equation is an enhancement factor of 9/4 due to the difference in the corresponding Casimir eigen- 
values. 

The numerical results for the adjoint representation are shown in the Fig. ^ On the left panel, 
we compare the fermion wave functions, (p^). when we use the modified BC (red dashed line) 
and CP vertices (black continuous line). 

On the right panel, we show the fermion dynamical mass ^adjl/'^)- We see that the infrared 
saturation of ^adj (p^) occurs for higher values compared to the values of ^{p^) in the fundamen- 
tal. More specifically, when the modified BC vertex is employed, one obtains ^adj(O) = 750 MeV, 
while for the CP vertex .-^adj (0) = 962 MeV. Clearly, due to the nonlinear nature of the gap equa- 
tion, the results found in the adjoint representation can not be reproduced from the fundamental 
solutions through a simple multiplication of the factor 9/4. For example, the value for ^adj (0) are 
clearly higher than 9/4^(0). 
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Figure 7: Left panel: The fermion wave function A^^jj {p^) in the adjoint representation. Right panel: The 
dynamical fermion mass ^adjl/?^)- 
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5. Conclusions 

We have shown that the quark gap equation can give rise to phenomenologically compatible 
results for the dynamical quark mass provided that (i) a complete non-AbeUan quark-gluon vertex 
(ghost-dressing function and quark-ghost scattering kernel) is introduced, and (ii) the recent lattice 
results for the gluon and ghost propagators are used. 

It is important to emphasize that the incorporation of the appropriate ingredients in the way 
dictated by the underlying symmetry, as captured by the Slavnov-Taylor identity for the quark- 
gluon vertex, furnishes to the kernel (especially in the intermediate region of the momenta) the 
required support, thus obviating the need to resort to additional (artificially induced) enhancements, 
of questionable field-theoretic origin. 
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